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Abstract. - We analyze the interplay between wetting and anchoring of nematic liquid crystals 
on disordering, e.g., rough substrates in the framework of the Landau-de Gennes theory, in 
situations of competing homeotropic and planar easy axes on the substrate and the nematic- 
isotropic (NI) interface. The phase diagram for azimuthally symmetric substrates is calculated. 
We identify two regimes - a strongly coupled regime, where the wetting transition coincides 
with an anchoring transition, and a weakly coupled regime, where the two are separated. The 
anchoring transition is first order and switches between homeotropic and planar anchoring. The 
two competing orientations are metastable over a broad parameter range. Hence such surfaces 
can be used to generate bistable surfaces. 



Nematic liquid crystals are fluids of aligned particles with a preferred orientation [1]. In the 
bulk, the direction of alignment is arbitrary, but surfaces break the isotropy of space and tend 
to align nematic fluids. This phenomenon, called surface anchoring, is of both fundamental 
and industrial interest, e.g., in liquid display technology [2,3]. From the practical point of 
view, it is useful to have surfaces which favor two orientations, such that one can effortlessly 
switch from one to another. Various strategies have been proposed to obtain such bistable 
surfaces [4-10]. Here, we analyze a mechanism that generates bistability on rough substrates. 



Experiments [11-13], theory [14], as well as simulations [15,16] have shown that rough 
or nanostructured substrates can reduce the order at the surface and even depress the NI 
transition in confined systems. If the disordering effect of a substrate is strong enough to 
nucleate a surface layer of strongly reduced order, the adjacent nematic fluid is oriented by 
that layer rather than by the bare substrate [12], and the direction of alignment may change. 
For example, Simoni et al. have observed experimentally that a porous polymeric substrate 
imprints a different orientation on a nematic film than a flat substrate made of the same 
material [17]. In their system, the alignment was homeotropic on the flat substrate, but planar 
on the porous substrate. This effect can be understood quite naturally if one assumes that 
the substrate favors homeotropic alignment, and the NI interface induces planar alignment. 

The question is whether the effect can be used to generate bistable surfaces. Rodriguez- 
Ponce et al. [18] have carried out density functional calculations for a model liquid crystal 
on a disordering substrate in a similar situation of competing, planar and homeotropic, an- 
choring axes. Encouragingly, they find a first order anchoring transition between planar and 
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homeotropic anchoring. The anchoring transition seems related to wetting, but the exact 
nature of the relation is not clear: The system also exhibits wetting transitions that are not 
connected to an anchoring transition, and even reentrant wetting. In the present work, we 
present an analysis of the problem in the general framework of the Landau-de Gennes theory. 
We clarify in detail the relation between the anchoring transition and the wetting transition, 
and identify different regimes of weak and strong coupling, which establish the context for 
the findings of Rodriguez et al. We find that the anchoring transition is first order, and 
that there exists a broad intermediate regime, where both anchoring orientations are at least 
metastable and the system can easily switch from one to the other. Hence these surfaces 
should be suitable candidates for bistable surfaces. 

The fact that competing easy axes may influence wetting phenomena has first been pointed 
out for nematic wetting layers by Sullivan and Lipowsky [19]. Director distortions in nematic 
wetting layer may lead to long-range interactions between the substrate and the NI interface. 
The consequences for the wetting phase behavior have been analyzed within the Landau- 
de Gennes theory by Braun et al. [20], using a formalism originally developed by Sen and 
Sullivan [21]. Within the same framework, Teixeira et al. [22] have studied a system with an 
isotropic wetting layer. This work is related to ours, but the system is much more complicated. 
Already the bare substrate free energy has two competing contributions: A term which is 
linear in the order parameter and favors surface order and homeotropic anchoring, and a 
quadratic disordering term which favors conical anchoring. As a result, a continuous surface 
driven anchoring transition from conical to homeotropic was found for a set of parameters 
with strong linear term, and a first order wetting driven anchoring transition from conical to 
planar for a second set with weak linear term. Here, we shall focus on the simpler case of a 
purely disordering substrate that favors homeotropic anchoring. 

Our starting point is the Landau-de Gennes theory of nematic liquid crystals, which is 
based on a free energy expansion in powers of a symmetric and traceless (3 x 3) order tensor 
field Q(r). The leading terms contributing to the bulk free energy are [1] 

F bulk = yrf 3 r{|Tr(Q 2 ) + |Tr(Q 3 ) + ^Tr(Q 2 ) 2 + ^d l Q jk d l Q jk + ^-d l Q lj d k Q k ^. (1) 

For simplicity, we neglect the biaxiality and approximate the order tensor by [23] Qij(r) = 
^S(r)(3rii(r)nj(r) — Sij), where S(r) is the local nematic order parameter, and n(r) the 
director, a vector of length unity describing the local direction of alignment. This is justified 
by the fact that homeotropically orienting surfaces do not induce biaxiality, and that the 
biaxiality induced by the NI interface is small [24,25]. Furthermore, we introduce the "natural" 
units S = -2B/9C, £ = 2^/(Li + L 2 /6)/3C S' 1 , and e = (3C/16) • S 4 £ 3 . The quantity S is 
the value of the order parameter in the nematic phase at coexistence, £ is the minimum width 
of a planar NI interface at coexistence, and e/£ 2 the corresponding interfacial tension. In the 
following, all order parameters, lengths, and energies, shall be rescaled by these units. This 
leaves us with two dimensionless parameters, 

t= \ A ^--l, and * L * (2) 

4 e 2 (Li +L2/6) V ; 

The parameter t is proportional to the distance to NI coexistence in the phase diagram - i.e., 
the temperature distance (t oc (T — Tjvj)) in thermotropic liquid crystals, or the chemical 
potential distance (t cx (/1 — ^ni)/Tni) m lyotropic liquid crystals. The parameter a charac- 
terizes the anchoring strength of a planar NI interface. At a > 0, the interface favors parallel, 
planar alignment, and at a < 0, it favors perpendicular, homeotropic alignment. 
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The resulting rescaled bulk free energy takes the form 

Fbuik = 3 J d 3 r {f + 9l + 52 } with f = S 2 ((S-l) 2 + t), (3) 

gi = ((V-S) 2 +a(n- VS) 2 ) + 4aS'((V ■n)(n-V5) + ^(nx V x nXVS)), 
.92 = S 2 ((3 + 2a)(Vn) 2 + (3 - a)(n ■ V x n) 2 I (3 + 2a)(n x V x n) 2 ) . 

The first term, f(S), describes the free energy density of a homogeneous system, the middle 
term accounts for the effect of order parameter variations, and the last term corresponds to 
the Frank elastic energy of a nematic phase with spatially varying director [1]. This term, 
allows to relate the parameter a to the experimentally accessible Frank elastic constants Ki 
of a material. For example, for MBBA, one has K 3 /K 2 = (3 + 2a)/(3 — a) « 3, i.e., a w 1.2. 

Next, we must determine the appropriate surface free energy. An isotropic surface intro- 
duces only one symmetry breaking vector, the surface normal no. This vector can be combined 
with Q to construct the surface energy terms that are compatible with the symmetry of the 
system [26]. The linear order term in Q, noQno = |5(nno) 2 — hS, favors nematic order, 
i.e., a surface free energy containing such a term will always be minimized by a nonzero value 
of S. Hence this term must vanish close to a truly disordering surface. To quadratic order in 
Q, one obtains three terms, from which one can construct the general expression 

F smi = J d 2 r / surf with / surf = WS 2 {\+l3n\+ 1 n^ and n\ := 1 - (nn„) 2 . (4) 

At W > and (3 > 0, the surface favors 5 = (disorder) and rz.ii = (homeotropic alignment). 
The parameter W > measures the disordering effect of the substrate, i.e., it's roughness, 
and the parameter (3 > characterizes it's orienting strength. For 7 < — (3/2, the surface 
has an additional preference for planar anchoring. In the following, we shall assume 7 = 
for simplicity. It is worth noting that the form ^ of / sur f already implies that the wetting 
transition must be second order, for symmetry reasons, unless it is coupled with an anchoring 
transition [27]. 

The total free energy is given by F = Fbulk + ^surf) with -Fbuik and F sur { given by Eqs. © 
and (HJ). Our task is to minimize this functional with respect to the profiles S(r) and n(r). We 
take the surface to lie in the (xy) plane, hence we can assume that the profiles vary only in the 
z direction. After parametrizing the director nasn = (n\\ cos0, n\\ sin0, n z ) with n 2 +n 2 z = 1, 
one checks easily that F as a function of n z and <j> is minimized by d<j>/dz = 0: The director 
does not vary in the azimuthal direction. Hence we are left with two profiles, S(z) and n z (z), 
which have to be determined such that they minimize the total free energy. 

We will now sketch a method that allows to solve this and similar problems very efficiently. 
We divide the order parameter S(z) profile in piecewise monotonic parts. (In our case, the 
whole profile was monotonically increasing). For each part, we rewrite the director profile as a 
function of the order parameter S, n z (z) = n z (S). After introducing new variables 9 = ln(5) 
and ip = 2 arcsin(n z ) for convenience, the bulk free energy ([3]) per surface area A can be 
written 

3^{/(e') + (|) 2 e 2 ^ 2 (^,^)} (5) 
with <I> 2 (V>, V') = 1 + ~(1 - cost/;) + ~ip' simp + 3 ^^ i)' 2 . 
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Fig. 1 - Examples of effective anchoring potentials E as a function of the anchoring angle 9 between 
the director and the surface normal. The parameters are W = 2.5, a = 1.2, (3 = 1., and t as indicated. 



It is first minimized with respect to q(z) for given ip(q). The Euler-Lagrange equations yield 
the integration constant f(e q ) — (dq/dz) 2 e 2q & 2 = const. = /oo, where /oo = /(Sncmatic) is 
the free energy density in the homogeneous nematic bulk. This can be used to derive an 
expression for dq/dz, which can be inserted into the total free energy, yielding 



F/A 



Vf^/A + 6 / dq e q v7(e«) 
J qo 



surf 



(6) 



with /gurf = We 2qo (l + /3cos 2 (?/>o/2) + -/ cos 4 (ip / 2)). The index oo stands for the bulk, the 
index for the surface, and V is the total volume of the system. The free energy ([6]) can 
now be minimized with respect to ip(q). A variational treatment yields the Euler-Lagrange 
equation 

^ ^77o = — - irrzrfl> - (1 + 7; — Hf(e q ) - /oo) ) — , (7) 



dip' 2 dip dip dip' 



2dq 



'I dip 1 ' 



<9$ 

with boundary conditions — — =0, = ln(5 ncma ti c ), 

dip' 9oo 



9$ 



and ee^V/Ceg)-/^ — 



9/ ( 



surf 



dip 



6e«V/(e 9o )-/oo^o,^) = 



9/, 



surf 



(8) 
(.9) 



Eq. ([U) (left) was obtained by minimizing (J6J) for arbitrary upper integration limit g max < , 
and then taking the limit q lrLa , x — > (/oo- One easily checks that both ip = and ip = tt are 
solutions of Eqns. ©-([I]). To calculate the anchoring potential for arbitrary anchoring angle 
6, we fix ipoo — tt — 29, calculate ip'^ from Eq. (|8|), and perform a straightforward integration 
of Eq. ([7]), starting at q^ and stopping as soon as Eq. © (right) is fulfilled. The resulting 
total free energy per area F/A = S gives the anchoring potential. At the extrema of S, the 
remaining boundary condition ([9]) (left) is automatically fulfilled. 

Some examples of anchoring potentials £(#) are shown in Fig. [T] As a rule, the anchoring 
potential always assumed it's minimum either at 9 = (homeotropic anchoring) or 9 = n/2 
(planar anchoring) [28] . Far from the NI coexistence and for weakly disordering substrates (low 
W), the orienting force of the substrate dominates and the effective anchoring is homeotropic. 
Close to the coexistence and for strongly disordering substrates (high W) , the main orienting 
force stems from the fluid layer with strongly varying order parameter close to the surface. 
In that case, the effective anchoring is planar. Fig. [2] shows two examples of phase diagrams 
in the (W,t) plane. At the transition lines, the surface is truly bistable, both homeotropic 
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Fig. 2 - Anchoring phase diagrams in the (W, t) plane (distance from NI coexistence (— t) vs. surface 
roughness W), for two values of the substrate orienting strength /3, at a — 1.2. The solid lines 
mark the phase transition, and the hatched areas indicate the parameter regions where both planar 
and homeotropic anchoring are at least metastable. The cross indicates the position of the wetting 
transition at t = 0. If one increases W at t — for small f3 (J3 < j3* = 1.16), one first encounters 
a first order homeotropic/planar transition, and then a second order wetting transition. For large 
(J3 > /?*), the homeotropic/planar transition triggers the wetting transition, which is then first order. 



and planar alignment are equally favorable. For practical purposes, it will often be sufficient 
if a state is metastable, i.e., if it corresponds to a minimum of F(9). Fig. [2] shows that the 
regions where both states are metastable (hatched areas) are quite large. Moreover, the energy 
barriers between the two states are small (~ 0.1e/£ 2 , see Fig. [J), hence switching between the 
two states is easy. 

For a more thorough understanding of this transition, we must relate it to the wetting 
transition. Thus we consider the two competing states, 9 = ir/2 and 9 — at NI coexistence 
(t = 0). The resulting free energy as a function of the surface order parameter, So, is 

f (l-S*o) 2 (l + 25 ) + (l + /3)PFS* 2 planar 

\ VTT^(1- S ) 2 (l + 2S Q ) + WS^ homeotropic. 1 ' 

If £ takes it's minimum at Sq = 0, the surface is wetted by the isotropic phase; otherwise, it is 
nonwet. For fixed anchoring angle, the wetting transition hence takes place at W* = 3/(1 +/3) 
in the planar case, and at Wj* = i^Jl + a in the homeotropic case, and it is continuous (critical 
wetting). In addition, the system may switch from the homeotropic state to the planar state. 
As expected, the wet surface always favors planar anchoring. 

The resulting phase behavior depends on the orienting strength (3 of the surface. If \f3\ is 
larger than a critical value |/?*|, the anchoring is homeotropic for all nonwet surfaces. The 
transition from homeotropic to planar anchoring coincides with the wetting transition, and 
is first order. For smaller the (first order) homeotropic-planar transition preempts the 
wetting transition, which is then continuous. In this case, the system switches to planar 
anchoring at a stage where the NI interface is not yet fully developed. Fig. [3] shows the 
wetting phase diagram for a — 1.2 and, more generally, the value f3* which separates the two 
regimes as a function of a. 

Rough, bistable, surfaces can be combined with conventional smooth or rubbed surfaces 
in a nematic liquid crystal cell device. As an example, we briefly discuss the phase diagrams 
of a nematic fluid confined between a bistable surface, and strongly orienting surfaces with 
fixed anchoring angle 9s — (homeotropic) or 9s — 7r/2 (planar). Since S = Soo in the 
cell, the only relevant contribution to the bulk free energy, ([3]), is the Frank elastic energy 



6 



EUROPHYSICS LETTERS 



Homeotropic 




Nonwet 




P" 


, Planar, Wet 

\ 

\ 

\ 




\ 

\ 

V 

Planar N 


^^-^ Nonwet \^ 



(b) 



1 2 

w 



Fig. 3 - (a) Wetting phase diagram at coexistence (t — 0) for a = 1.2. Solid lines correspond to first 
order, dashed line to continuous transitions. For |/3| < |/3*|, the planar/homeotropic phase transition 
and the wetting transition are separated; for \/3\ > |/3*|, they coincide, (b) Critical value (3* where the 
wetting and the anchoring transitions meet as a function of a. Here, results are also shown for the 
case a, (3 < 0, where the surface favors planar alignment and the NI interface aligns homeotropically. 



52, and we can write the total free energy per area of the system as F cc \\/A = X(#r) + 3(3 + 
2a) 5^,(^5 — Or) 2 /D. Here Or is the anchoring angle on the rough surface, and D is the 
thickness of the film. This must be minimized with respect to Or. The phase diagrams for a 
system with the same surface parameters as in Fig. Q] are shown in Fig. 2] The system can 
assume two different configurations - a homogeneous configuration with fixed director angle, 
= Or throughout the system, and a hybrid aligned nematic (HAN) configuration where the 
director slowly rotates from the planar to the homeotropic orientation in the z-direction. The 
transition lines between the different states can be calculated to a very good approximation, 
if one assumes that Or can only take the values Or = or Or — ir/2 (dashed lines in Fig. |4j. 
Hence the rough surface really acts like a two-state surface. 

To summarize, we have analyzed the interplay of wetting and anchoring on disordering 
substrates within the Landau-de Gennes theory in situations of competing, homeotropic and 
planar, anchoring axes at the substrate and the NI interface. We have identified different 



(a) 





HAN 


Planar 


8 R =it/2 



(b) 



1000 HAN 



Homeotropic 



e„=o 



0.2 0.4 0.6 0.8 1 
-t 



0.1 
-t 



Fig. 4 - Phase diagrams of a nematic liquid crystal confined between a bistable surface and a strongly 
anchoring surface with anchoring angle 6r — n/2 (a) and 9r = (b). HAN denotes a "hybrid aligned 
nematic" state where the director rotates from homeotropic to planar alignment. The two other 
states correspond to homogeneous planar or homogeneous homeotropic alignment. The thick solid 
lines indicate the exact transition lines, dashed lines an approximative result as explained in the text. 
The parameters of the bistable surface are W — 2.5, a — 1.2, (3 = 1. 



F. Schmid and D. Cheung: Bistable anchoring 



7 



substrate regimes - one "strongly orienting" regime where the anchoring transition coincides 
with the wetting transition, and one "weakly orienting" regime where it preempts the wetting 
transition. The anchoring transition is first order, but the wetting transition must be contin- 
uous, if it is not coupled with an anchoring transition. Finally, we have discussed how this 
effect can be used to design bistable surfaces, which favor two distinctly different orientations, 
and given an example how such a surface could be integrated in a nematic liquid crystal cell. 
The present study is a mean field study, the effect of fluctuations has been disregarded. Two 
types of fluctuations will renormalize the surface potential: Director fluctuations become in- 
creasingly important close to the NI transition [29], and fluctuations of the NI interface [30] 
become important at complete wetting. It will be interesting to study the influence of these 
effects in future work. 

This work was funded by the DFG and the EPSRC. 
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